Monte Carlo study of the scaling of universal correlation lengths in three-dimensional 

0(n) spin models 
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Using an elaborate set of simulational tools and statistically optimized methods of data analysis 
we investigate the scaling behavior of the correlation lengths of three-dimensional classical O(n) spin 
models. Considering three-dimensional slabs S 1 x S 1 x R, the results over a wide range of n indicate 
the validity of special scaling relations involving universal amplitude ratios that are analogous to 
results of conformal field theory for two-dimensional systems. A striking mismatch of the n — > oo 
extrapolation of these simulations against analytical calculations is traced back to a breakdown of 
the identification of this limit with the spherical model. 
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I. INTRODUCTION 

The concept of scaling, the observation that singu- 
lar observables vary in a scale-free manner according to 
power laws when the driving parameter of a transition 
(temperature, magnetic field, ...) is tuned towards a criti- 
cal point, has since the first observations beea a key ingre- 
dient of the theory of critical phenomena™! Exploiting 
the symmetry of scale-invariance, forming the geometri- 
cal basis for the power-law behavior in the vicinity of a 
critical point, through the idea of real-space renormal- 
ization, scaling theory can be mapped on the behavior 
of finite systems near the transition point of the bulk 
system in the limit of diverging system sizes, the-,|t|hcr- 
modynamic limit. This finite-size scaling (FSS)tirclEl oc- 
curs with scaling exponents generically linked to the ex- 
ponents that govern scaling in the bulk system. Thus 
the apparent weakness of finite system size that hampers 
simulational approaches actually turns out to be their 
intrinsic strength, when exploring FSS means exploring 
thermal scalingElEl. 

The significance of scaling theory for the understand- 
ing of critical phenomena becomes quite exposed in 
the context of conformal field theory (CFT) for two- 
dimensional systemsEi. In the course of exploiting the 
additional invariances of conformal symmetry one is able 
to split the critical point partition function of a lattice 
system into a sum over contributions from all the scaling 
variables present in a specific model. Consider a critical 
system on a L x L' lattice with toroidal boundarsj-eondi- 
tions; then the partition function decomposes asBE3: 
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where c is the central charge of the considered theory, / 
the bulk free energy per unit volume, 5 — L' /L, A = LL', 
and the sum runs over the whole content of scaling op- 
erators with dimensions x n . Thus, the knowledge of the 
operator content of a theory in connection with the cor- 
responding scaling dimensions is equivalent to an "exact" 
solution of the model on finite lattices. 



Two Dimensions - A particular example of a scaling 
relation in two dimensions that can be derived assuming 
conformal invariance of critical point entities concerns 
the two-point function in the limit of L' — > oo. It it gen- 
erally sufficient to assume translational, rotational, di- 
latational_and invcrsional invariance to imply conformal 
invarianceO; homogeneity, isotropy and scale invariance 
alone suffice to uniquely fix the critical, connected two- 
point function of an operator </> in the infinite plane up 
to an overall normalization factor: 



(4>(zi, z\)(j)(z 2 , z 2 )) c = (zi - z 2 ) x (zi - z 2 y 



(2) 



where z\ , z 2 are complex co-ordinates parametrizing the 
plane. Then, one uses the logarithmic map 
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to wrap the complex plane around an infinite length 
cylinder S 1 x M of circumference L with co-ordinates 
w = u + iv, where v measures the polar angle along 
S 1 and u the longitudinal direction along R. Assum- 
ing conformally covariant transformation behavior of the 
(primary) operator 0, one arrives at^n expression for the 
two-point function on the cylindertd: 



(4>(w 1 ,W 1 )(j)(w2,W2)) 
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2tt\ ( n , 27T. , „ 2tt . , 

— I I 2 cosh — (m - u 2 ) ~ 2 cos— [v\ - v 2 ) ). (4) 



In the limit of large longitudinal distances |ui — u 2 \ 3> L 
and v\ = v 2 , one is left with a purely exponential drop 
with a correlation length 
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Thus, utilization of conformal invariance yields a finite- 
size scaling relation including the amplitude, which is in 
contrast to renormalization group theory that usually 
gives the scaling exponents and only certain amplitude 
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ratios, but not the amplitudes themselves. Since this 
result emerges from a field-theoretic description of sta- 
tistical mechanics that does not take into account the 
microscopical details of the system, it is expected to be 
universa&3. Note, however, that this proposed universal- 
ity goes beyond the usual notion of an universal quantity 
and comprises three different aspects: (i) the correlation 
length of a given operator should be the same within the 
associated universality class of models; (ii) when looking 
at different operators, on the other hand, the form of Eq. 
(||) should be left unchanged, all operator-dependent in- 
formation being condensed in the scaling dimension x; 
(iii) finally, even when looking at models of different uni- 
versality classes, all that should change are the scaling 
dimensions (and the definition of </>) , the validity of Eq. 
(||) being untouched. Property (i) implies t le "hyper- 
universality" relation of Privman and FisheiO. In the 
following, we will refer to the whole extent of aspects (i)- 
(iii) exceeding the usual notion of universality with the 
term "hyperuniversal" . A corollary that is of importance 
for transfer matrix calculations that use an unnormalized 
(quantum) Hamiltonian results from taking the ratio of 
the correlation lengths of two primary operators, for ex- 
ample the densities of magnetization and energy which 
are usually primary for spin models: 
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Because of the independence from the overall amplitude 
1/27T of Eq. (||) this relation might still stay valid when 
changing the geometry in a way such that only this over- 
all amplitude is altered. In terms of universality this con- 
stitutes a weaker form of the aspect (i) above, namely 
universality of amplitude ratios instead of amplitudes 
themselves; we will refer to this weaker property as (i') 
in the following. 

A suitable test-bed for these general field-theory re- 
sults is, of course, given by the exactly solvable two- 
dimensional Ising model. Using Eq. (||) and the generic 
relations between scaling dimensions and the conven- 
tional critical exponents: 







(7) 



giving x a = 1/8 and x t = 1 for the two-dimensional 
Ising model, one arrives at a ratio x e /x a = 8. A di- 
rect evaluation of the spin-spin correlation length in the 
Onsager-Kaufman framework gives, as the leading term 
in the scaling series, £= = iLL/7r = L/(2ir^), in agreement 
with the CFT result00ll3. The same holds true for the 
leading scaliua amplitude of the energy-energy correla- 
tion functiontSl, £ e = L/2n. Both amplitudes have also 
been evaluated nuiacrically to high precision in a Monte 
Carlo (MC) studyll3, resulting in perfect agreement with 
Eq. (|). 

A possible alteration of the S 1 x K situation, namely 
changing the boundary conditions along the ^-direction 



from periodic to antiperiodic has also been treated within 
the CFT framework, exploiting the fact that in the case 
of the ferromagnetic nearest-neighbor Ising model the 
antiperiodic boundary corresponds to the insertion of 
a seam of an£iferromagnetic_honds along this boundary 
line. This calculation yieldsEl'EJ: 
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again in good agreement with Monte Carlo dataEl. Note, 
however, that this last relation, in contrast to Eq. (JsJ) , is 
specific to the Ising model and the special choice of the 
densities of magnetization and energy as operators and 
thus is not "hyperuniversal" in the sense of properties 
(ii) and (iii) presented above. 

The amplitude-exponent relation Eq. (|^) for two- 
dimensional systems has been checked analytically or 
numerically and found valid for an impressive series of 
further rpedels like the Eeitts model and its percola- 
tion lipajill3, the XY modeled, the syppetric eight-vertex 
modellS, and quantum spin modelsEj to name only the 
most prominent. 

Three Dimensions - On leaving the domain of two- 
dimensional systems towards higher dimensions, the 
wealth of exact field theoretic calculations is instantly 
reduced to severe scarcity. The conformal group co- 
incides with the set of holomorphic functions in the 
special case of spatial dimension d = 2 and is thus 
infinite-dimensional as a group. For d > 3, unfortu- 
nately, it reduces to a simple Lie group with dimension 
D < (d + l)(d + 2)/2 for any Riemannian, connected 
manifold. As a consequence, only in two dimensions the 
postulate of conformal invariance is restrictive enough 
for a classification of the operator contents of the differ- 
ent universality classes and thus an exact solution of the 
critical theories within the limits of field-theory assump- 
tions. For d > 3, on the other hand, the implications of 
the finite-dimensional conformal-group symmetry reach 
hardly beyond the consequences of plain renormalization 
group theory exploiting dilatational invariance. However, 
since inversional symmetry is still present, a transforma- 
tion like Eq. (^) stays conformal in higher dimensions, 
now connecting the spaces R d and S^ 1 x R. Applied 
to the two-point function one arrives at a scaling rela- 
tion analogous to Eq. (||), namely £y = R/x, cp. Rcf.Ej, 
which contains the d = 2 result as a special case assuming 
L = 2wR, R being the radius of S 1 . Since primarily 
of operators is a priori not well defined for d > 3, it is, 
however, unclear for which operators this relation should 
hold. A numerical analysis for this geometry, which has 
to cope with the fact that S' d ~ 1 for d > 3 is a truly curved 
space and thus hard to regularize by discrete lattices, will 
be presented in a separate publicationEj. 

On the other hand, the toroidal geometry S 1 x . . . x 
S 1 x R, which is much more convenient for numerical 
simulations, is not conformally flat and thus no CFT pre- 
dictions exist for this case. In spite of this theoretically 
unfavorable situation a transfer matrix calculation for the 
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Hamiltonian limit of the three-dimensional Ising r 



on the geometry S 1 x S 1 



T 2 x R by HenkeJ^ 



rendered results still comparable to the situation for the 
x R geometry. For the ratios of leading scaling 
amplitudes of correlation lengths for different boundary 
conditions (be) he found 



£ CT /£ e = 3.62(7) for periodic be, 
£ CT /£ e = 2.76(4) for antiperiodic be. 



(9) 



A comparison with the (inverse) ratio of the correspond- 
ing scaling dimensions, 



(1 - a)/v _ 2(i/d- 1) 
vd 



2.7264(13), (10) 



[3/v vd — 7 

(cp. Table | and Eq. (Q)) showed that even though 
the original expectation to possibly find agreement in 
the case of periodic boundary conditions as in the two- 
dimensional case was not met, the data are consistent 
with the relation Eq. (^) for the unorthodox case of an- 
tiperiodic boundary conditions. Note that one has to 
compare ratios in this case, because the quantum Hamil- 
tonian used in the calculation is defined only up to on 
overall normalization constant. This result is in qual- 
itative .a&reement with a Metropolis MC simulation by 
WestorJ 2 ^, who found ratios £<x/£e of about 3.7 for peri- 
odic and 2.6 for antiperiodic boundary conditions, re- 
spectively. Considering these striking observations it 
seems interesting to check whether this behavior is just 
a coincidence or special feature of the Ising model or in- 
stead indicates a general property of critical models on 
this special three-dimensional geometry. 

The rest of the paper is organized as follows. In Sec. II 
we introduce the general class of models we want to ex- 
amine and present the way we are going to discretize the 
three-dimensional geometry T 2 x R. We discuss simula- 
tion methods, observables, estimators for measurements 
and parameters of the simulations. In Sec. Ill we outline 
the statistical tools used for the data analysis. It is quite 
hard to extract high-precision information about corre- 
lation lengths from MC simulation data; we will thus 
discuss the special path of data analysis we are going to 
proceed along and present details of the statistical tools 
used there for. This tool-set is "calibrated" with simu- 
lations of the two-dimensional Ising model, where exact 
results for comparison are available. In Sec. IV we dis- 
cuss the results for the correlation lengths ratios of our 
simulations for the Ising, XY and (generalized) Heisen- 
berg_rnodels. Our results, already briefly announced in 
Ref.EJ, confirm Henkel's findings on a high level of accu- 
racy. Furthermore this behavior seems to carry through 
for the whole class of O(n) spin models and is thus far 
from being a "numerical accident" . In Sec. V we try to 
rank our numerical findings in the context of the clas- 
sification of universality presented above. The type of 
the model considered enters not only via a variation of 
the scaling dimensions, but also influences the overall 
prefactor of Eq. (||). Sec. VI is devoted to the discus- 
sion of the relation of our finite-n results to the spherical 



TABLE I: Literature estimates for the critical exponents v 
and 7 of the three-dimensional Ising model. 

Method 
g-expansio: 
e-expansionl 



series 
series 
MC^, 
MC^ 

M£.'t 
MCB 



12 



weighted mean 



0.6304(13) 


1.2396(13) 


0.6305(25) 


1.2380(50) 


0.6315(8) 


1.2388(10) 


0.63002(23) 


1.2371(4) 


0.6289(8) 


1.239(7) 


0.6301(8) 


1.237(2) 


0.6303(6) 


1.2372(13) 


0.6298(5) 


1.2365(5) 


0.63005(18) 


1.23717(28) 



model, which is connected to the limit n — ► oo of the class 
of 0(n) spin models. The classic identification of both 
models seems to break down as soon as (multi-point) 
correlation functions are considered. The final Sec. VII 
contains our conclusions. 



II. MODELS AND SIMULATION 

Throughout this paper we consider classical, ferromag- 
netic, zero-field, nearest- neighbor, O(n) symmetric spin 
models with Hamiltonian 



H = 



(11) 



The underlying lattice is taken to be simple cubic with 
dimensions L x x L y x L z . Special cases of this class 
of models include the Ising (n = 1), XY (n — 2), and 
Heisenberg (n = 3) models. This Hamiltonian has the 
advantage of representing a whole class of models with 
critical points in three dimensions, tuned by the single 
parameter n. According to the T 2 x R geometry we set 
L x = L y and apply periodic or antiperiodic boundary 
conditions in the x and y directions. In both cases we use 
periodic boundary conditions in the ^-direction to elim- 
inate surface effects that are also absent in the L z — > oo 
case assumed in Eq. (||). To reduce effects of finite size 
in z-direction one has to ensure that L z 3> £|| , a concrete 
rule will be given below. 

In view of the problem of critical slowing dawn, we 
use the Wolff single cluster update algorithmic for all 
O(n) model simulations, cp.EJ. The adaption of this 
update procedure to the case of antiperiodic boundary 
conditions along the torus directions is straightforward 
if one exploits the above mentioned equivalence of an 
antiperiodic boundary to the insertion of a seam of an- 
tiferromagnetic bonds along the boundary line for the 
case of nearest-neighbor interactions. Considering the 
Ising model c*_ alternatively, embedded Ising spins for 
n > 1 modelsE3, this means: adjacent spins interacting 
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FIG. 1: Typical spin configurations for the two-dimensional Ising model on strips of size 20 x 382. (a) periodic boundary 
conditions; (b) antiperiodic boundary conditions. Note that the visible geometric clusters differ from the stochastic clusters of 
the cluster update algorithm. 



antifcrromagnetically are connected with a bond obey- 
ing the Swendsen-Wang probability p = 1 — cxp(— 2(3 J) 
in case of opposite orientation and are left unbonded in 
case of identical orientation. This rule exactly reflects 
the change in energy compared to the ferromagnetic case 
and thus trivially satisfies detailed balance. 

The main observables of our simulations are the con- 
nected correlation functions of the densities of magneti- 
zation and energy: 

G=(xi,x 2 ) = (o-(xi) • <r(x 2 )> - {a) ■ (tr>, , , 
^( Xl ,x 2 ) = (e(x 1 ) e (x 2 ))-(e)(e). 

We define the energy density as a local sum over the 
nearest neighborhood x' of a spin x (x' nn x) : 

<^ = - J 2 £ (13) 

x' nn x 

the factor 1/2 ensuring that E = ^ x e(x). It is straight- 
forward to construct a bias-reduced estimator for the case 
of (x 2 — Xi) || e 2 , corresponding to the correlation length 
£ = £h: first, taking advantage of the translation in- 
variance of the systems along the z-axis established by 
a periodic boundary, one can average over the "layers" 
i = \z2 — z\ | = const, r-^fo improve on that consider a 
"zero-mode projection"EJ, i.e. define layered variables 

W = tV £ (14) 

* x' ,z—z 

where Ot = c f or et denotes the times series of MC mea- 
surements, and consider the estimator 

T 

= fY,T z £ O t ( Zl )O t (z 2 ) 

-(^EEaw)'. < 15 > 

where T denotes the length of the MC time series. This 
estimator obviously does not directly measure G C: ", but 
inspecting the continuum form Eq. reveals that the 
deviation stemming from transversal cross-correlations 



entering the estimator declines exponentially with in- 
creasing longitudinal distance i and thus becomes irrele- 
vant for the long-distance behavior we are interested in. 
Numerical investigations confirm that these considera- 
tions stay correct when passing to three dimensionsta. In 
the large-distance regime zero-mode projection reduces 
the variance of correlation function estimates by a factor 
inversely proportional to the layer volume L x L y . Note 
that the given estimator for the disconnected part (O) 2 
has a bias that vanishes as l/T in the large-T limit. 

As mentioned above, periodic boundary conditions in 
z-direction eliminate surface effects associated with this 
direction, but still effects of finite L z will trigger devia- 
tions from the L z — > oo limit assumed in Eq. (|^). Inspect- 
ing the form of Eq. (^) in the limit of distances i £n 
one expects longitudinal correlations according to: 

G C: " (i) oc e-^n + e-^'-^o , (16) 

i.e. the exponential decay is superimposed by an expo- 
nentially increasing part. Thus, using too small values of 
L z results in an effective underestimation of correlation 
lengths. In order to satisfy L z 3> £n in a systematic way, 
i.e. to keep this effect away from the region of clear signal 
for measuring the correlation lengths, and assuming lin- 
ear scaling of correlation lengths according to £|| = AL X , 
one has to keep the ratio L z /^ = L Z /AL X fixed and 
therefore has to scale L z proportionally to L x . Simu- 
lations for the case of the two-dimensional Ising model 
show that these finite-size effects are negligible compared 
to the statistical errors for L z /£\\ > 10 and lengths of 
time series of about 10 6 to 10 7 measurementstj. Adding 
a safety margin the longitudinal system sizes for the 
simulations in three dimensions where chosen such that 
~ 15, the scaling amplitude A being estimated 
from a simulation of an "oversized" system. Since £ CT > £ e 
for all models under consideration, the amplitude A a of 
the spin-spin correlation length scaling is significant for 
the satisfaction of this condition. Note that from Eq. (|l5| ) 
increasing L z also has the positive side effect of improv- 
ing the statistics of the correlation function estimation. 

In order to judge the efficiency of the used cluster up- 
date algorithm and to ensure reasonable usage of com- 
puter time we evaluated integrated autocorrelation times 
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Ti n t using a binning technique!^]. The strong asymme- 
try of the model lattices reduces the average size of 
clusters and thus Wolff's cluster update algorithm does 
not perform as good as on (hyper-)cubic lattices, result- 
ing in increased autocorrelation times. Since measure- 
ments of G c '" are computationally expensive compared 
to update steps, but the statistical gain vanishes with 
increasing r; n t, measurements were done with frequen- 
cies of about 1/rjnt- Approaching the low-temperature 
phase, antiperiodic boundary conditions in the torus di- 
rections produce a spatially stable boundary of the geo- 
metric clusters along the antiferromagnetic seam, which 
in turn enforces a second boundary along the z direc- 
tion. This results in a further reduction of the average 
cluster size compared to the periodic boundary case. Fig. 
|l| shows typical configurations for the case of the (two- 
dimensional) Ising model. 

III. DATA ANALYSIS 

Having sampled correlation functions according to 
Eq. ( |l5|) and assuming the functional form G C: "(i) = 
a exp (— ) + b, we refrain from using instrinsically un- 
stable non-linear three-parameter fits and resort to the 
following estimator instead, 



£o(i) = A 



In 



-i -l 



(17) 



which eliminates the additive and multiplicative con- 
stants a and b above. Note that it is not allowed to 
assume b — a priori for time series of finite length, cp. 
Ref.EJ. Apart from stability considerations this approach 
allows for computational simplifications, since correlation 
functions can be sampled irrespective of normalization 
and the biased estimation of the disconnected part (O) 2 
can be dropped. In addition, Eq. ( |l7| ) simplifies the dis- 
tinction of the long-distance part of the correlation func- 
tion from the short-distance region: as the explicit two- 
dimensional expression Eq. (||) implies, exponential decay 
will only occur asymptotically, but with deviations de- 
caying themselves exponentially; apart from that, lattice 
artefacts that are not reflected in the continuum expres- 
sion Eq. (El) additionally distort the short-distance behav- 
ior. Fig. g shows an example plot of the spin-spin cor- 
relation length estimates £o-(«) for the Ising model. The 
transition from the short-distance region that should not 
be used for the final estimate to the purely exponential 
long-distance behavior is clearly visible. The parameter 
A in Eq. ( |l7| ) can be used to tune the signal-noise ratio 
for the correlation length estimate; increasing A dramat- 
ically reduces the apparent statistical fluctuations in £(i), 
cp. Fig. ||. Note, however, that the reduction of variances 
for individual distances i is accompanied by an increase 
of cross-correlations between estimates for adjacent esti- 
mates, so that the error of an average over a region of 
distances becomes minimal for a value A clearly below 
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FIG. 2: Correlation length estimates according to Eq. (JTTj) 
and O = cr for a 30 2 x 382 Ising system with periodic bound- 
ary conditions for two choices of the typical distance A. The 
plateau regimes collapse if both ordinates are scaled identi- 
cally. 



its allowed maximum. As a compromise, we use A 
for both estimators ^(i) and £ e (i). 
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Naive estimates for the statistical errors (variances) of 
complex, non-linear combinations of observable measure- 
ments like the estimator Eq. ([Tt]) are extremely biased 
due to two effects: even for quite sparse measurements 
with frequencies around 1 /T- mt successive elements of the 
time series are still correlated, generically leading to sys- 
tematic underestimation of variances. This effect is be- 
ing eliminated by the grouping together o£jpeasurements 
to sub-averages of length fi ( "binning" )EJ, which leads 
to an asymptotically uncorrelated time series of length 
T" = T I fi used in the further process of error estimation. 
For the production-run time-series the bin size was cho- 
sen to regularly include several thousand measurements, 
which is far in the asymptotical regime. Secondly, the 
strong non-linearity of estimators like Eq. ( |l7j ) forbids 
the use of the usual formula for the standard deviation 
of a set of measurements. A common solution to this 
problem is the use of the Gaussian error propagation for- 
mula, which, however, only uses a lowest order Taylor 
series approximation to the functions and assumes Gaus- 
sian distribution of the mean values, i.e. long enough time 
series for all observables. A far more general ansatz is. 
given by resampling techniques such as the "jackknife"c3 
that apply to a quite general set of probability distribu- 
tions and capture function non-linearities exactly. The 
jackknife variance and bias estimators mimic the brute 
force error estimation method of comparing k completely 
independent MC time series of lengths T' and applying 
the naive estimates: removing single elements (i.e. bins) 
of a single time series of length T 1 one by one results in 
T" time series of length T' — 1, e.g. for the correlation 
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function estimates: 



(18) 



resulting in jackknife-block estimates for the correlation 
length of: 



in 



Gf w (i)-G(.)(*-A) 



n -1 



G {s) (i + A)-G s (i) 

few = ^ £*«(*•). 



(19) 



Then the jackknife estimate of variance is given by: 



vAR(&)) = E (e w (») - 4)« 



(20) 



Note the missing factor of 1/ (T" — l) 2 as compared to the 
usual variance estimate which accounts for the trivial cor- 
relation between the T' jackknife-block estimates. One 
can show that this estimator, apart from the reweighting 
prefactor (T" — 1)/T', is strictly conservative, i.e, devia- 
tions from the true variance are always positiveo. Simi- 
larly, the resampling scheme provides an estimate for the 
bias of estimators, namely: 



BIAS(£(i)) = (T'-l)(£ (0 (t)-|(i)), 



(21) 



and thus offers a bias corrected correlation length esti- 
mate as = T'i(i) - (T' - Since in non- 
pathological cases the bias of an estimator vanishes with 
increasing length of the time series, the jackknife bias 
estimate provides a good check for whether the consid- 
ered series are long enough to neglect bias. A jackknife 
error estimate for these bias-corrected estimators is pos- 
sible iterating the jackknife resampling scheme to second 
ordered. 

Since Eq. (0) gives a vector of estimators for the cor- 
relation length instead of only a single one, an improved 
final estimate can be achieved by an average over the 
£(«). However, as for example Fig. |^ reveals, only a cer- 
tain range of distances i — i m i n , ■ ■ ■ , W is suited for 
this purpose, where the lower bound i m ; n results mainly 
from small-distance deviations as reflected by Eq. (|^), 
whereas the large distance bound i max cuts off the re- 
gion where the signal of exponential fall-off drops be- 
low the size of statistical fluctuations, so that error es- 
timates become inaccurate and eventually the estimator 
Eq. (|lj]) becomes maldefined due to negative arguments 
of the logarithm. Conventionally, averaging over the es- 
timates for i = i m i n , . . . , ?max would be done with 
weights Ui oc 1/ a 1 (£(«)) that minimize the theoretical 
variance of the mean value. This prescription, however, 
neglects correlations between the individual estimates. 
Note that cross-correlations between adjacent estimates 



are quite large, not only because large scale fluctua- 
tions of the correlation functions arc dominant, but also 
since the used estimator Eq. ([l7]) explicitly introduces 
such correlations increasing in range with increasing A. 
As a simple variational calculation shows, for the case 
of correlated variables to be averaged over, one has to 
choose the weights according as 



Oik 



(22) 



in order to minimize the variance of the mean value. 



Here, T e 



■>, P 



1, denotes the co- 



variance matrix of the £(i). T itself can be estimated 
within the jackknife resampling scheme as: 



CORRy = CORR(£(i),£(j)) 



T' — 1 



'j 1 E (£« w -4) w) (lw (j) - 4) ui) ■ (23) 



T 



The fact that, considering Eq. p2|), variance and co- 
variance estimates directly influence the final results for 
the correlation lengths, gave the motivation for the quite 
careful statistical treatment presented above. 

Finally, the selection of the regime i = i m in, • ■ ■ , imax 
can, besides the obvious eyeball method, also be done in 
a way based on statistical criteria. Interpreting the aver- 
age over the as a fit of the estimated values to 
the trivial function /(£) = £ = const, the systematic de- 
viations from the plateau regime for very small and very 
large distances i should be clearly reflected in quality-of- 
fit parameters. Thus, looking at the x -distribution, 



x 2 = 



= E ( 24 ) 



will be a good criterion for judging the "flatness" of 
the plateau regime i m i n , ■ ■ • , 'max included in the average. 
Again, as an estimator for the covariance matrix one 

can use the jackknife expression CORRy . Then finding 
the optimal region of distances for the average is equiv- 
alent to the optimization problem \x 2 /g — 1| — » min, 
with g — i maK — i 1U in = p — 1 denoting the number of 
degrees of freedom of the fit. However, this ansatz of 
optimization bears some uncertainties: minimizing the 
distance of x 2 1 9 from 1 supposes that the optimal choice 
includes estimates whose dispersion around £ is ex- 
actly reflected by the estimated variances. In view of the 
jackknife's tendency to overestimate errors it might be 
more favorable to minimize \x /g\ itself. Furthermore, 
considering the statistical nature of the data, the abso- 
lute minimum of \x 2 1 9 ~ 1| or IxVffl sometimes happens 
to be an isolated fluctuation, far apart from the bulk of 
next-to-optimal solutions. Finally, this optimization pro- 
cedure tends to result in minimal values for very small 
regime sizes p since the fit becomes trivial for very small 
numbers of points; this, however, conflicts with another 
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FIG. 4: Finite-size scaling of the energy-energy correlation 
length of the three-dimensional Ising model with antiperiodic 
boundary conditions. The other scaling plots look similar; we 
show the worst case. The fit was done to the functional form 
Eq. (§1). 



0.1583(17), in excellent agreement with the exact results 
A a = 4/tt ps 1.27324 and A t = l/2vr « 0.15915, cp. Eq. 
(§). For the case of antiperiodic boundary conditions 
we arrive at A a = 0.42410(30) and A t = 0.07984(38), 
compared to CFT results of A a = 4/3vr w 0.42441 and 
A c = 1/4tt w 0.07958, cp. Eq. (§). 



FIG. 3: Sections of X 2 /3(*min, imax) for the spin correla- 
tion length of an Ising system, (a) {% 2 /g \ imin = 25}; (b) 
{xVtflw = 110}. The ' 'wavy" structure results from 
A = 4 in Eq. ©. 



possible goal of optimization, namely the minimization 
of the overall variance of the final result. To circum- 
vent these problems we resort to considering the whole 
two-dimensional distribution x 2 /s(*min, W)- It is char- 
acterized by a rather flat plateau regime for intermediate 
values of i mul and i max and steep increases at the bound- 
aries, cp. Fig. ||. A good recipe for the determination of 
bounds is then given by first choosing a preliminary i m - ln 
well above the steep ascent for small i; then a plot like 
Fig. ||(a) allows to determine the upper bound z max - Fi- 
nally, a plot of {x 2 1 g | i max = const} determines the final 
lower bound z m in, cp. Fig. ||(b). 

To test the methods of data analysis described in this 
section we performed simulations of the two-dimensional 
Ising model. Using a series of systems with L x = 
5, . . . , 20 and finite-size scaling fits including an effective 
higher-order correction term of the form £(£ x ) = AL X + 
BL^, we find for the leading correlation lengths scal- 
ing amplitudes A a u final estimates for the case of peri- 
odic boundary conditions of A a — 1.27374(81) and A e = 



IV. RESULTS: AMPLITUDE RATIOS 



Let us now turn to the three-dimensional geometry 
T 2 x R and the determination of amplitude ratios ac- 
cording to Eq. (||) . We report the results of simulations 
for the O(n) spin models for n = 1, 2, 3, and 10. 

Ising Model - Simulations of the Ising model were 
done at an inverse temperature given by a high-precision 
MC estimate of the bulk critical coupling in three 
dimensions^, /3 C = 0.2216544(3). We use a tempera- 
ture reweighting technique to check for .the .influence of 
the uncertainty of (3 C on the final resultsoa. We find it 
completely negligible compared to the statistical errors 
for the case of the Ising model. To enable a proper FSS 
analysis including sub-leading terms we performed simu- 
lations for transverse system sizes L x = 4, 5, . . . , 20, 25, 
and 30, scaling L z accordingly. Adapting the frequency 
of measurements to the autocorrelation times, about 
2 x 10 6 and 8 x 10 6 nearly independent measurements 
were recorded for the systems with periodic and with 
antiperiodic boundary conditions, respectively. Collect- 
ing the final estimates £ for the correlation lengths one 
ends up with a scaling plot like that shown in Fig. |^. 
The scaling behavior is quite linear, however, as plots 
of the amplitudes £/L x reveal, corrections to the purely 
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TABLE II: Literature estimates for the inverse critical tem- 
perature /3 C and the critical exponents v and 7 of the three- 
dimensional XY (n = 2) and Heisenberg (n = 3) models. 



resultE 



n Method 



7 



e-expansior 

32 

series; 

serie; 
serie; 
serie 
MC^ 
MCl! 
MCl 
MCi 
MCO 



0.45419(3) 

0.45406(5) 
0.45420(6) 
0.4542(1) 
0.454165(4) 

0.45421(8) 
0.454170(7) 



0.6680(35) 

0.677(3) 

0.67166(55) 



0.670(2) 
0.672(1) 
0.6723(11) 



1.3110(70) 

1.327(4) 

1.3179(11) 



1.319(2) 
1.316(3) 
1.3190(22) 



weighted mean 0.454167(3) 0.67179(42) 1.31839(82) 



e-expansio 



MCi! 

MCi; 

MCEI 



0.6929(1) 

0.69305(4) 

0.6929(1) 

0.693035(37) 

0.6930(1) 

0.693002(12) 



0.7045(55) 

0.712(10) 

0.716(2) 

0.706(9) 

0.7036(23) 

0.704(6) 

0.7128(14) 



1.3820(90) 

1.400(10) 

1.406(3) 

1.390(23) 

1.3896(70) 

1.389(14) 

1.399(2) 



weighted mean 0.69301(1) 0.71129(98) 1.3998(16) 



linear scaling behavior are clearly resolvable, cp. Fig. |^. 
As an aside, Fig. ||(b) additionally shows jackknife bias 
corrected estimators according to Eq. (|2l|); for the given 
length of time series bias effects of our estimator Eq. ( |l7| ) 
can clearly be neglected. 

Returning to the two-dimensional case for a moment, it 
is easy to see the source for the leading correction term 
in the correlation length scaling. In the framework of 
conformal field theory the effect of lattice discretization 
as well as the influence of non-linearity of scaling fields 
that increase with the distance from criticality (i.e. the 
thermodynamic limit in our case) can be included, in con- 
siderations using conformal perturbation theorylliJ. A for- 
mal perturbation expression for the spin-spin correlation 
length including the effect of a perturbing operator cou- 
pled with strength ctk is to first order given by 



2tt 



27rafc(Cii,i — C fco) 



(25) 



where the perturbing operator has dimension Xk and the 
coefficients C n kn result from the operator product ex- 
pansion (OPE). One findso that to lowest order the only 
non- vanishing amplitude belongs to an operator that cor- 
responds to the breaking of rotational symmetry by the 
square lattice as compared to the continuum solution. 
It has dimension Xk = 4 leading to 1/L corrections, 
in agreement with the first-order expansion of the exact 



2tt 

T 



2tt 



1 /2tt 
768tt I ~L 



(26) 



A similar effect will be present in the three-dimensional 
systems, but the correction exponent can no longer be 
evaluated analytically. Fig. || shows that the sign of 
the leading correction term is unchanged in three dimen- 
sions for the systems with periodic boundary conditions, 
whereas it is reversed for the systems with antiperiodic 
boundary. This stays true for the other 0(n) spin models 
discussed below. To account for corrections to scaling we 
fit the correlation lengths data to the functional form 



£(L X ) = AL X + BL* 



(27) 



treating the correction exponent re as an additional fit 
parameter. Due to the presence of higher-order correc- 
tions, however, the resulting values of re have to be taken 
as effective exponents, that will in general differ from 
Wegner's correction exponent u). Therefore we decided 
to keep re as a parameter, despite of existing field-theory 
estimates for u>, cp.El. We take into account the effect of 
neglecting higher-order correction terms by successively 
dropping points from the small L x end while monitoring 
the quality-of-fit parameters x 2 /.9 or Q to find a com- 
promise between fit stability and precision of the final 
amplitudes A. The range of sizes L x used is indicated by 
the range of the solid lines in Fig. ||. Our results for the 
scaling amplitudes and their ratios as listed in Table ET3 
and the ratio of scaling dimensions according to Eq. (|10|) 
show precise agreement in the sense of Eq. (||) for the 
case of antiperiodic boundary conditions and clear mis- 
match for a periodic boundary. Thisis. in agreement with 
the results of HenkelE3 and WestonEa, but at a level of 
accuracy that makes a casual coincidence very unlikely. 

XY Model - The XY model is, as well as the Heisen- 
berg models, accessible to cluster update methods using 
the embedded cluster representatiorH, which we made 
use of. The simulations were performed at the coupling 
j3 c = 0.454167(3), which is an average of recent literature 
estimates, cp. Table H Using the same transverse system 
sizes L x — L y as for the Ising model, but adjusting the 
lengths L z according to the different correlation length 
amplitudes, we took between 1 x 10 6 and 16 x 10 6 mea- 
surements, using measurement frequencies around l/ri nt 
as above. Fig. ^ shows the amplitude scaling plot of the 
spin-spin correlation length for periodic boundary condi- 
tions. The additional curves are results of a temperature 
reweighting analysis, trying to judge the effect of critical 
coupling uncertainties. The precision of the data is well 
illustrated by the fact that, reweighting our results to the 
minimum and maximum estimated critical couplings, re- 
spectively, cited in Table O, results in a variation of the 
scaling curves far beyond the range covered by the re- 
maining statistical errors. Nevertheless, reweighting to 
the lcr-range inverse temperatures (3 C — A/? and [3 C + A/3 
as given above triggers deviations at most comparable 
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FIG. 5: Scaling of the amplitudes £ a / e /L x for the Ising model. The solid lines show fits to the function Eq. (J27|) ; (a) and 
(b) show correlation lengths for the systems with periodic boundary conditions, (c) and (d) for the case of an antiperiodic 
boundary; (b) additionally contains bias corrected estimates according to Eq. (Ell). 



to the error estimates of the statistical analysis. The 
intermediate maximum of the curve for /3 m i n , however, 
might be an artefact indicating that /3 m in is already too 
far away from the simulation temperature to allow for 
reliable reweighting. The effect of temperature variation 
is generally observed to be smaller for the antiperiodic 
boundary systems; furthermore, it is more important for 
the case of the spin-spin correlation length since here 
statistical errors are clearly smaller than for the energy- 
energy correlation length estimates. Thus, Fig. ^ shows 
the largest effect observed. 

Fitting the final correlation length results £ CT / e to the 
functional form Eq. (p7j), we arrive at the final estimates 



for the leading amplitudes given in Table [II . Comparing 
these to the ratio of scaling dimensions resulting from the 
averaged critical exponent estimates of Table [lj and Eq. 
(|l0|), we again find Eq. (^) confirmed for antiperiodic 
boundary conditions only; this behavior is obviously not 
specific to the Ising model. 



Heisenberg Model - The n = 3 Heisenberg model case 
is treated analogously to the XY model. Table || gives 
the critical parameter estimates used for the simulations 
and comparison. With statistics similar to that for the 
n = 1 and n = 2 cases, the simulations confirm the find- 



ings of the Ising and XY models, cp. Table III for details 



For the case of the energy-energy correlation length of the 
systems with periodic boundary conditions the gathered 
statistics did not suffice for a stable non-linear fit includ- 
ing corrections according to Eq. fl27| ) . We thus performed 
a simple linear fit dropping the correction term. This, 
however, results in an error estimate which is not quite 
realistic and, furthermore, induces a systematic under- 
estimation of the amplitude since one expects B e < 0, 
cp. Fig. ||(b). From the results of the other models this 
effect is estimated to be about 2o~-3o~ in magnitude. 

O(10) Model - To gain additional evidence and in 
order to facilitate considerations concerning the n — > oo 
limit, giving a clear picture of systematic dependencies 
on the parameter n, we also simulated the n = 10 gen- 
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FIG. 6: Amplitude scaling of the spin-spin correlation length 
of the XY model with periodic boundary conditions. The 
spread curves show results of temperature reweighting for j3 c — 
A/3 = 0.454164, [3 C + Af3 = 0.454170, /3 min = 0.45406 and 
max = 0.45421. 



TABLE III: FSS amplitudes of the correlation lengths of O(n) 



spin models on the T 



geometry. 



Model 




JJcIlULllC DC cLLlLlLJcllOLllL, UL 




A 

A a 


U.oloo^OZ ) U.ZoOy4^oU j 


Ising 


At 


U.2zoz(lo J U.UoOOl(ol ) 


A 1 A 
Aa j rl e 


O.yjyJyjyOVJ } Zi. (OU^lOj 




"T E j <"G 


2.7264(13) 




A 

A a 


7^400^0^ 941 1 1(^7} 


XY 


\ 

A e 


U.loyy^lO J U.UoZo^lo j 


A 1 A 
Aa j s±c 


3 Q71 1^9^ 9 Q^ni'47^ 




/ Xa 






A 

A a 


n 7oc\f<9if'^A\ n 9/i/i^9(' p ii ~\ 




A e 


0.16966(36) 0.0793(20) 


Heisenberg 


A a /A e 


4.2478(92) 3.085(78) 






3.0891(79) 






0.671107(59) 0.25865(46) 


n = 10 


A, 


0.1350(23) 0.07096(107) 


A a /A f _ 


4.971(83) 3.645(55) 






3.615(70) 



eralizcd Heisenberg model. Since, of course, in the past 
much less effort has gone into the investigation of the 
0(n) model with n > 4, there are quite few estimates 
of the critical coupling. We thus here use a single high- 
temperature series estimate of (3 C — 2.42792(8)t3. The 
implementation of the Wolff cluster update algorithm has 
to cope with the technical intricacy of generating pseudo- 
random numbers equally distributed on a hyper-sphere, 
see Appendix A for details. Due to this complication 
we only simulated systems up to a transversal size of 
L x = 20 and reduced the number of measurements to 
2 x 10 6 . The critical exponents for comparisntL-given by 
a plain average over some recent estimates^EHtEj, are: 



v = 0.8713(75), 7 = 1.721(14). 



(28) 



Table III shows again agreement between amplitude and 
exponent ratios only for the case of antiperiodic bound- 
aries. Note that, as critical exponent estimates become 
rare with increasing n, the correlation length ratio esti- 
mate already reaches the precision of the scaling dimen- 
sion ratio estimate. Checking the influence of the critical 
coupling uncertainty we find it only important compared 
to statistical errors in the case of the spin-spin correla- 
tion length for periodic boundary systems; the results 
reweighted to 0± = /3 C ± A/3 are A~ = 0.670805(56) and 
A+ = 0.671432(65), respectively. This, however, does 
not noticeably influence the error of the ratio estimate, 
since here the error of the estimate of A e , which is much 
larger, is dominant. 

We thus find the linear amplitude-exponent relation 
Eq. (JgJ) confirmed for several spin models in three dimen- 
sions with the peculiarity that one has to insert a seam 
of antiferromagnetic bonds along the redirections to re- 
store the two-dimensional situation. 



V. RESULTS: "META" AMPLITUDES 

Comparing our results for the three-dimensional ge- 
ometry S 1 x S 1 x M to the CFT conjecture for the case 
of two dimensions, we are interested in the respective 
ranges of validity in terms of the classification of uni- 
versality aspects given above in the Introduction. The 
fact that our simulations of the isotropic lattice repre- 
sentation of the O(n) universality classes give results in 
agreement with the strongly anisotropic quantum Hamil- 
tonian representation used by Henkel in his trans, 
trix calculations for the case of the Ising modell 
indicates that the considered amplitude ratios— are uni- 
versal, i.e. (i') holds. Apart from that, HenkeEa explic- 
itly checked for universality of amplitude ratios by the 
insertion of an irrelevant perturbing operator and found 
it confirmed for both cases of boundary conditions. How- 
ever, strictly speaking, there is no proof of universality 
for the cases n > 1. The universality aspect (i) above, 
i.e. universality of the amplitudes themselves, could not 
be checked in Henkcl's calculations, because the quantum 
Hamiltonian is only defined-up to an overall normaliza- 
tion constant. Yurishche\CJl23 considered the behavior of 
an anisotropic Ising model and found varying correlation 
lengths amplitudes on variation of the ratios of couplings 
in the different directions. This, however, is no argument 
against amplitude universality since anisotropy is repre- 
sented by marginal instead of irrelevant operators. On 
the other hand, amplitude ratios stay universal even with 
respect to those marginal perturbations, in consistency 
with Henkel's strongly anisotropic Hamiltonian limit cal- 
culations. In fact it has been argued that for all systems 
below their upper critical dimension correlation length 
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scaling amplitudes are universal quantities^. 

Having found very good agreement in three dimensions 
between ratios of correlation lengths and scaling dimen- 
sions according to Eq. (|^) for the case of antiperiodic 
boundary conditions, it is interesting to see what the 
overall, operator- independent, "meta" amplitude A ac- 
cording to: 



Ccr/e — ^ff/e-^i — 



A_ 

Ccr/e 



(29) 



that was A = 1/2tt for two-dimensional periodic systems, 
cp. Eq. (||), becomes in three dimensions, in particular 
whether it is again model independent. Since our results 
for the spin-spin correlation lengths are always more pre- 
cise than those for energy-energy correlation lengths, we 
use ^o- to determine A. The estimates for the spin-spin 
scaling dimension x a resulting from the corresponding 
estimates of bulk critical exponents v and 7 listed in Ta- 
bles I and H and Eq. (§|) are x a = 0.5182(4) (Ising), 
x a = 0.5188(9) (XY), x„ = 0.5160(17) (Heisenberg) , and 
x a = 0.512(12) (n = 10), respectively. Thus, inserting 
our results for A a listed in Table III, we obtain for the 
"meta" amplitudes A(n): 



A AfjX(j — 



0.12278(43) Ising 

0.12510(37) XY 

0.12622(49) Heisenberg 

0.1325(30) n = 10 



(30) 



These values can additionally be compared with an an- 
alytical result that is available for the case of the spherical 
model, which is commonly believed tp-be identical to the 
n — > 00 limit of the O(n) spin modeltJ. Agaia-using the 
Hamiltonian formulation, Henkel and WestorO'Ea found 
that the amplitude exponent relation Eq. (^) is valid for 
the spherical model on S 1 x S 1 x 1 for both kinds of 
boundary conditions, periodic and antiperiodic. This is 
due to the fact that the quantum Hamiltonian factorizes 
into a set of uncoupled harmonic oscillators. The ampli- 
tude A for the case of antiperi odiq boundary conditions 
was found to be A ~ 0.13624E3l3. Plotting this value 
together with the finite-n results of Eq. ( p0| ) shows an 
apparently smooth variation of the "meta" amplitudes 
with the order parameter dimension n, the eyeball ex- 
trapolation of the finite-n values to 1/n — > matching 
the spherical model result, cp. Fig. 0(a). Facing this vari- 
ation, the hypothesis of a "hyperuniversal" amplitude 
A(n) = A that does not depend on n, as was the case for 
the two-dimensional systems, can be clearly ruled out. 
Thus, type (iii) universality of the classification above 
gets broken when passing from two to three dimensions. 
The matching of the finite-n values with the universal 
spherical model amplitude, on the other hand, indicates 
universality also of the finite-n amplitudes and thus uni- 
versality of type (i) above. 

Even without a scaling law of the type Eq. (0) be- 
ing valid for the case of periodic boundary conditions, 
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FIG. 7: (a) "Meta" amplitudes A for antiperiodic boundary 
conditions according to Eq. ( po| ) as a function of the order 
parameter dimension n; (b) The same combination A a x a for 
periodic boundary conditions according to Eq. (|3l|). 



one can nevertheless plot the corresponding combination 
A a x a for this case also, as is illustrated in Fig. 0(b). The 
values are: 



AfjXn 



0.4240(17) Ising 

0.3912(7) XY 

0.3719(12) Heisenberg 

0.3439(78) n = 10 



The corresponding value for the spherical model is given 
by A a x a ~ 0.3307, cp.EaE3. The finite-71 values again run 
smoothly into the spherical model limit. 



VI. THE LIMIT OF INFINITE SPIN 
DIMENSIONALITY 

While the finite-n amplitudes of Fig. fit well to the 
spherical model result, this is not the case for the corre- 
lation lengths ratios themselves. From inspection of Fig. 
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FIG. 8: Correlation lengths ratios as function of the order pa- 
rameter dimension n for periodic and antiperiodic boundary 
conditions. 

|§| the smooth variation of correlation length ratios for fi- 
nite n does not fit at-all to the spherical model result of 
Henkel and WestontOS that gives a ratio A a /A e = 2 for 
both, periodic and antiperiodic boundary conditions. By 
eyeball extrapolation one would instead expect the am- 
plitude ratios to reach values around 4 for antiperiodic 
and around 5^ for periodic boundary conditions in the 
limit n — > oo. And indeed, accepting the validity of a 
linear amplitude-exponent relation according to Eq. (|^) 
for the case of antiperiodic boundary conditions and us- 
ing the usual relations for the connection between scaling 
dimensions and bulk critical exponents, namely Eq. (Q), 
one would expect x a — 1/2 and x e = 2 since (3 = 1/2, 
v = 1 and a — — 1 for the spherical model. The resulting 
ratio x c /x a = 4 perfectly agrees with the eyeball extrap- 
olation of our finite- n data. However, by inspection of the 
energy-energy correlation function in the Hapajltonian 
limit and using factorization arguments, HenkelEJ conjec- 
tured x e = 1 instead, resulting in the ratio A a jA t = 2, in 
contrast to the relation Eq. (R) . Taking into account the 
obvious agreement of eyeball extrapolation and spherical 
model calculation for the amplitudes A(n) that were cal- 
culated from the spin-spin correlation length amplitude 
as A{n) = AfjXfj , cp. Fig. |7|, it becomes obvious that the 
mismatch is entirely due to the behavior of the energy- 
energy correlations. Note also that, since the specific heat 
is constant in the low-temperature phase of the spherical 
model in three dimensions, interpreting this as an effec- 
tively vanishing specific-heat exponent a' = leads to 
an effective energetic scaling dimension x' f = 1. This, in 
fact, implies a violation of the scaling relation Eq. (]?]), 
which is of the hyper-scaling type, for the case of the 
spherical model. 

Puzzled by this striking mismatch, we performed a 
roughly explorative MC simulation directly in the spher- 
ical model, which rendered results in qualitative agree- 
ment with an amplitude ratio of A a /A e = 2 as suggested 



by the analytical calculation. Then, it is natural to ask 
whether there is a contradiction with Stanley's result on 
the equivalence of the p_— > oo limit of the 0(n) model 
and the spherical modelO, which has been, after some da- 
bate over mathematical subtletiesO, rigorously provenO. 
The precise statement that can be proven is the identity 
of the partition functions or, equivalently, free energies of 
the two models in the thermodynamic limit for the whole 
temperature range, even independent of the order of tak- 
ing the limits n — > oo and N — > oo (the thermodynamic 
limit). Since multi-point correlation functions do not fol- 
low from the (source-free) partition function, this does 
not say anything about the behavior of these functions 
in those two models. A direct calculation in the spherical 
model, cf. Appendix B, results in a simple factorization 
property of the long distance behavior of the connected 
energy-energy correlation function for all temperatures 
in one and two dimensions and in the high-temperature 
phase down to T c in three dimensions. If the four-point 
function of the spherical-model spins is denoted by Cijki, 
one has: 

Cii+l j j+l ~ = Ci jCi+1 j+1 + Cj j + lCi+1 j 

-^2Cfj, |j-i|-oo, (32) 

where Cy are the corresponding two-point functions. 
This confirjaas Henkel's results for the Hamiltonian 
formulationEJ on more general grounds. 

Considering the n — > oo limit of the 0(n) model, on 
the other hand, reveals that the connected part of the 
energy-energy correlation function vanishes in the first- 
order saddle-point approximation that is being used for 
the comparison of the two models, cf. Appendix C. This 
is in agreement with general considerations for the large 
n model by BrezinEl For the case of the one-dimensional 
spin chain, the connected energy-energy correlation func- 
tion even vanishes exactly for all n, so that one can rule 
out an agreement of the two limits to higher order of 
the steepest-descent expansion in this case. Thus the 
mismatch of finite-n extrapolations and spherical model 
results of Fig. || has some well-defined mathematical rea- 
son. 

Starting from the observation that the curves of Fig. 
|§| for the amplitude ratios seem to be quite parallel as a 
function of (finite) n for the both kinds of boundary con- 
ditions, we also plotted the collapsed ratio ^."j^' that 
should be unity if the amplitude-exponent relation Eq. 
(||) holds true. Inspecting Fig. || this is, according to our 
above results, of course the case for antiperiodic bound- 
ary conditions. Moreover, and a priori somewhat unex- 
pected, this ratio seems to be also quite constant for the 
case of a periodic boundary, stabilizing around a value 
compatible with 4/3 within statistical errors. Note that 
the exceptionally small error of the value for n — 3 (the 
Heisenberg model) and its apparent deviation towards a 
larger ratio is due to the impossibility to fit the n = 3 
energy-energy correlation lengths to a scaling law includ- 
ing a correction term as mentioned in Sec. IV. Statisti- 
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FIG. 9: Matching of correlation lengths ratios A a /A e and 
inverse scaling dimension ratios x t /x a for the two kinds of 
boundary conditions as a function of the order parameter di- 
mension n. The horizontal lines show fits to a constant as 
discussed in the text. 



cally, the data are consistent with a fit to a constant 
(Q = 0.08), and perfectly so when dropping the n = 3 
point (Q = 0.4). 

In view of this observation one might argue that the 
asymptotic scaling relation Eq. (^3j) in three dimensions 
has to be replaced by a generalized ansatz of the form 



R- 



(33) 



with an overall, model independent factor R that de- 
pends only on the boundary conditions and happens to 
be just 1 for the case of an antiperiodic boundary. For the 
amplitude scaling law this would lead to an asymptotic 
form 



^ /e (n)=R^-L a 



(34) 



cp. Eq. (|^). Accepting such a generalized ansatz, a least- 
squares fit of the collapsed ratios of Fig. |to a constant 
R gives R = 1.0037(45) for antiperiodic boundary condi- 
tions, underlining the validity of the original amplitude- 
exponent relation Eq. (||), or alternatively Eq. (^3|) with 
R = 1, for this case. For the periodic-boundary systems, 
on the other hand, we arrive at R — 1.3546(76) (omitting 
the n = 3 point), indeed statistically consistent with the 
conjectured value of 4/3. 

This somewhat diminishes the at first sight apparently 
exceptional importance of choosing antiperiodic bound- 
ary conditions in three dimensions. Taking into account 
the smooth amplitude variation of Fig. |^(b) the same uni- 
versality statements hold for periodic and for antiperiodic 
boundary conditions. 



VII. CONCLUSIONS 

We performed extensive MC simulations for several 
representatives of the class of O(n) spin models. Con- 
centrating on the geometry of three-dimensional slabs 
S 1 x S 1 x 1 we found a simple inversely linear relation 
between the leading scaling amplitudes of the correlation 
lengths of the magnetization and energy densities and the 
corresponding scaling dimensions valid to high accuracy 
for the Ising (n = 1), XY (n = 2), Heisenberg (n = 3), 
and n = 10 generalized Heisenberg models for antiperi- 
odic boundary conditions along the torus directions. This 
is the analogue of the CFT result in two dimensions with 
periodic boundary conditions applied. There is evidence 
for the universality not only of amplitude ratios (type 
(i') of our classification in the Introduction), but also of 
scaling amplitudes themselves (type (i)). To definitely 
decide the question whether universality in the sense (ii) 
above, i.e. condensation of all operator dependent in- 
formation in the scaling dimensions, is present, further 
operators would have to be considered. Independence, 
apart from changes in the scaling dimension, of the scal- 
ing amplitudes from the model under consideration, i.e. 
type (iii) universality, is explicitly broken for three di- 
mensions as compared to the two-dimensional case: we 
find a smooth variation of the overall "meta" amplitudes 
A(n) — Aa-(n)x a (n), depending on the order-parameter 
dimension n. It might be interesting to consider further 
classes of models, such as for example Potts models, to 
see whether any of these properties are specific to the 
O(n) spin model class. 

Considering the deviation of the periodic boundary 
correlation lengths ratios from the corresponding inverse 
scaling dimension ratios, the validity of a scaling law of 
the form Eq. (^|) can be definitely ruled out for this case. 
Generalizing this ansatz with an overall factor R depend- 
ing on boundary conditions as in Eq. (p3|), however, we 
find it fulfilled also for the case of periodic boundaries 
with a factor R independent from n and taking a value 
compatible with 4/3. In view of that, the fact that R = 1 
for the case of antiperiodic boundary conditions might 
be rather a coincidence than a "deep" physical property. 
Taking into account that in two dimensions the corre- 
sponding prefactors are specific to the operators consid- 
ered, cp. Eq. (|§|), makes it probable that a similar behav- 
ior occurs in three dimensions, destroying type (ii) uni- 
versality. It might be interesting to analyze the behavior 
of correlation lengths in the four-dimensional geometry 
S 1 x S 1 x S 1 x R to check whether a scaling law of the 
generalized form Eq. ( |33"| ) can be retained and if so, how 
the factor R depends on the dimensionality of the lattice. 

Trying to match our finite n results with analytical 
calculations for the spherical model we found a strik- 
ing mismatch of the data concerning energy-energy cor- 
relations. Inspecting the four-point functions directly in 
the spherical model and the 0(n — > oo) model limit we 
find that both results do not match to first order of the 
saddle-point approximation in general dimensions and to 
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all orders in one dimension. Thus, the idea of equiva- 
lence of the two models has to be limited to its original 
extent, namely the identity of partition functions in the 
thermodynamic limit. Quantities not directly related to 
the partition function, like multi-point correlation func- 
tions, do not necessarily have to coincide. Further work 
has to be done to possibly evaluate exactly the correla- 
tion lengths ratios in the n — > oo limit for both sorts of 
boundary conditions. 

Since, still, there is no explanation of the findings con- 
cerning the correlation lengths ratios for finite n in terms 
of a field-theoretic or otherwise exact approach, we would 
like to encourage further research in this direction. 
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APPENDIX A: EQUAL DISTRIBUTION OF 
RANDOM NUMBERS ON A HYPER-SPHERE 



Consider a probability density in polar co-ordinates 
f(4>,9) equally distributed on the 2-sphere S 2 , i.e.: 



m&) d^dg 

sin 9 d(/)d9 



= const. 



(Al) 



Factorizing /(</>, 9) = p((j>) q(9) = const • sin 9, and taking 
into account the normalization condition J dilf((j>, 9) — 
1, one finds: 



f(<P,e)=p(cf > )q(9) 



1 

2^ 



1 



(A2) 



Pseudo-random number generators usually generate 
numbers equally distributed in the unit interval [0,1]. 
How docs this transform to an arbitrary distribution? 
Let a random variable z be distributed with a density 
g(z) and transform according to z' = uj(z); the density 
h(z') then follows from the equation 



g(z) dz = h(z') dz' = h(uj(z)) oj'(z) dz. 



(A3) 



Thus, for random numbers z equally distributed in [0, 1] 
the transformation 9 = arccos(l — 2z) gives the de- 
sired distribution q(9) — ^sin9. This form is being 
used for the simulations of the n — 3 Heisenberg model. 
For general polar co-ordinates in R n , x\ ~ rcos(9i, 
X2 = r sin 9\ cos 92, up to x n — r sin#i • • • sin0 n _i, where 
< 9i < 7r, < #„_i < 2ir is understood, the volume 
element is given by: 

dV = r n - 1 sin"" 2 6»i sin"" 3 9 2 ■ ■ ■ sin6>„- 2 dr [J d9 l7 

(A4) 



so that one has for the factors f^'(0i) of an equally dis- 
tributed density /(0 X , . . . , 9 n . x ) = JJ, f {t) &Y 



7(71 — i — 1) 



sin"^ 1 ^, i<n-l, 



/(«- 1 )(0 n _ 1 ) = -, 

Z7T 



(A5) 



with normalization factors j(k) = 0Fr(^±i)/r(| + 1). 
Thus, for Zi equally distributed in [0, 1] the transforma- 
tions Zi(9i) are given by: 



7(n 



1) 



d9i sin" 



9i, (A6) 



for i < n — 1. The integrals can be evaluated analyt- 
ically for each 9i. There is, however, no closed form 
expression for the inverse transformation 9i(zi) that is 
needed to generate random vectors equally distributed 
on the hyper-sphere S n . The trivial workaround so- 
lution of sampling equally distributed in the hyper-cube 
L" = [—1, 1] x • • • x [—1,1], discarding the complement 
L n \B n and projecting the remaining points on the sphere 
S n , suffers from asymptotically vanishing efficiency, 
since the ratio of used to discarded volumes vanishes with 
increasing n exponentially as 7r"/ 2 /2"r(^- + 1). We thus 
resorted to a numerical inversion of Zi(9i) using interpo- 
lation between the pre-calculated points of a binary tree. 



APPENDIX B: ENERGY-ENERGY 
CORRELATION FUNCTION IN THE 
SPHERICAL MODEL 

Consider the spherical model of Berlin and Kac0 con- 
sisting of "spins" ei € K with the constraint: 



N 



(Bl) 



where TV denotes the number of lattice sites. For ease of 
reference we use the notation of the original paper here; 
thus, the 6i are not to be confused with the local energy 
densities defined above in Eq. ([l3j). The Hamiltonian is: 



H 



(ij) 



(B2) 



Using the Fourier representation of the (^-constraint Eq. 
(Bl) the partition function can be written as: 



Zn — 



A 

~2ni 



N 



ao+ioc 



d.S 



dei • • • dejy 



QO — ZOO 

' (ij) 



(B3) 



choosing ao such that the singularities in s of the in- 
tegrand are excluded from the integration volume. An 
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ensures the correct normalization of the integral mea- 
sure and K — j3J denotes the coupling. Diagonalizing 
the quadratic form with eigenvalues Xj via an 

orthogonal transformation = ^ , Vij yj , the Gaussian 
integration over the can be performed: 



fdyi--- dy N exp[- ^(s - KXj)y 

3 

[-5 EH 



Tr^exp 



ifA, 



(B4) 



so that, 

Z N = A N ^ N l 2 2Ke^ N ^ 2K ^- / dz 

1 N 1 

xexp[iV2^z--]rin(z--A J )], (B5) 

3=1 

where s = 2Kz. This expression can be evaluated in the 
saddle point limit N — > oo depending on the distribution 
of the eigenvalues Aj for a given lattice. Now consider 
the two-point function, 

dj = foe,-) = £ v *rV 3S (y r ys) = VirVjr , (B6) 



where the last equality foll ows from the symmetry of the 
partition functio n E q. (|B3| ). Compared to the Gaussian 
integration Eq. (B4) the insertion of a factor y 2 in the 
integrand gives an additional factor of: 



2(s-KX r ) AK{z-\X r ) 



(B7) 



The corresponding integral over ZrGan also be evaluated 
in the saddle point approximationO. Now, analogously, 
consider the four-point function: 

Ciju = (etCjekei) = VirVj a VktVi u (y r y s ytyu) ■ (B8) 



Here, again, only paired occurrences of the y m survive 
due to the inversion symmetry: 



Cijkl = Y V irVjrVkrVl r (y*) + ^ Vi r Vj r VkaVls (y. 



Ivl) 



r^s 



J2 VirV js V kr V ls (y 2 y 2 ) + £ V lr V JS V ks V lr (y 2 y 2 ). (B9) 



The insertion of yf under the Gaussian integral gives 
an additional factor of 3/[4(s - KX r ) 2 ] = 3/{16K 2 (z - 
X r /2) 2 }, whereas y 2 y 2 gives l/\16K 2 (z- X r /2)(z- X s /2)], 
so that the diagonal terms left out in Eq. ( |B9|) are rein- 
serted: 

Cijkl — J2 r _ s {VirVjrVksVls + VirVjsVkrVls 

+V ir V js V ks Vi r ) (y 2 y 2 ) (BIO) 



Now performing the z-integration of Eq. (B5) in the sad- 
dle point limit N — > oo is equivalent to just inserting the 
saddle point value z = z s for the factors given above, 
whenever a normal saddle point exists. As Berlin and 
Kac have shown, this is the case for all finite temper- 
atures in one and two dimensions and for T > T c in 
three dimensions; in the low-temperature phase, the sad- 
dle point "sticks" to its value for T = T c . Then, the 
four- point function simply factorizes, so that, comparing 
Eq. ( |B10| ) to the expression Eq. (B6) for the two-point 
function it is clear that: 



Cijki — CijCki + CikCji + CuCjk, 



(Bll) 



and, finally, considering the connected energy-energy cor- 
relation function, one has: 



a 



ii+l j j+1 



~ Cfi+i — Ci j C l+ i j+i + Ci j + \Ci+i , 



2C 2 V \j 



(B12) 



so that the energy-energy correlation function is triv- 
ially rela ted to the spin-spin correlation function. Note 
that Eq. ( Bll ) would follow from Wicks's Lemma for the 
Gaussian model. This especially confirms the factor-two 
relation x e /x a = 2 between the corresponding scaling 
dimensions derived by Henkel using transfer matricesEH. 
The factorization property can also be seen in the grand- 
canonical formulation of the spherical model, the "mean" 
spherical modelo, where the hard constraint Eq. (Bl) is 
being replaced by its thermodynamical average, so that 
one can leave out the problematic z-integration above. 
There has been some debate over the coincidence of the 
thermodynamic limit-of the two models, which is now 
believed to be scttlcdta. 



APPENDIX C: ENERGY-ENERGY 
CORRELATION FUNCTION IN THE LIMIT OF 
INFINITE SPIN DIMENSIONALITY 

The treatment of the partition function of the O(n) 
model in the n — > oo limit is quite analogous to that of 
the spherical model, cp.E3. For the comparison of the 
n — > oo limit with the spherical model the constraint 
<Tj • <Tj = 1 of Eq. (|Tl| ) has to be replaced by <Ji ■ cri = n. 
We write the partition function of the model as: 



Z%\K) = A 



(n)" 
N 



da 



(1) 



■•do- 



(n) 
N 



exp[x£5>^], 



(CI) 



(ij) " 



where A^ ensures the correct normalization. Rewriting 
the ^-constraints to the Fourier representation, one now 
has to introduce N variables {ti}, arriving at: 



Z%\K) = A%> {-- 



K 
2~i 



N 



-t-oc 



da\ 



(i) 



•da 



16 



+200 



dti • • • dtjy exp(Kn U 



seems only to hold to leading order of the saddle-point 
approximation. 



(C2) 



In the broken, low-temperature phase Eq. (C3) has to 
be replaced by 



u=l 



(ij) 



Interchanging the order of integrations one is again left 
with integrals of Gaussian type that are easily solved 
transforming the spin variables orthogonally according 



to a^' = J2j V%j yf ■ Note that the transformation is 
symmetric in the component index v of the spins. The 
calculation is given in more detail for the case of a one- 
dimensional chain below. Here, we again consider the re- 
lation between two-point and four-point correlation func- 
tions. We take the two-point function to be 



(C3) 



where the last equation for any v = 1, . . . ,n follows from 
the 0(n) symmetry of the model in the unbroken, high- 
temperature phase. Using the same arguments of Gaus- 
sian integration as for the case of the spherical model, 
the four-point function: 



1 



Cijki = — 5-((<Ti • er j)(cfc • (Ti)), 



(C4) 



again decomposes in terms of the diagonal variables y. 
as: 

-1 V V ri V rj V tk V tl (yM 2 yl u) 
n z * — 4 



(") 



Cijki 



\Y,V H V sj V sk V r i (yM 2 y 



(C5) 



In the saddle point limit, which now corresponds to 
n — > oo, this expression factorizes in terms of two-point 
functions as: 

Cijki = CijCki H — CikCji H — CuCjk, (C6) 
n n 

so that the "mixed" terms are suppressed with 1/n. This 
asymmetry results from the preset pairing of the spin 
component indices [i and v in the four-point function. As 
a consequence, the connected part of the energy-energy 
correlation function: 

Cji+l jj+l ~ = — C'i jCi+l 3+1 H Cjj + lCji+l: 

(C7) 

vanishes in the first-order saddle-point approximation. 
Thus, any non- vanishing contributions that are to be ex- 
pected from our numerical results, have to come from 
sub-leading terms in the steepest-descent expansion. The 
correspondence of the n — > oo limit to the spherical model 



Cij — — ((Ti ■ a 

n 



(C8) 



so that the factorization property of the four-point func- 



tion Eq. (C6) becomes 



iikl S ^ii^kl H <^ih ^ 



n J n J 



(C9) 



and again the connected part of the energy-energy cor- 
relation function is 0(l/n), vanishing in the first-order 
saddle-point limit. 

For the case of an one-dimensional lattice the first- 
order saddle-point approximation is exact as can be 
checked by explicit calculation. Consider an open chain 
of 0(n) spinsL3. The pa rtit ion function is given by the 
general expression Eq. (CI) with the nearest-neighbor 
sum fi • Cj replaced by the one-dimensional ex- 

pression J2i ° 'i ' 0i+i < Following Stanley^, we factor out 
the integration over the last spin er/v, which has the form: 

z m(k) = jl y . . .y da w . . . da (n) j +io ° du 

xexp[u^(n-^a (l/)2 )] exp[K ^ c v a {v \ (CIO) 



where c„ 



'jv-i- 



Inserting the unity factor 



exp[Kao(n— J2 V a ^ )] an d choosing a sufficiently large 
to exclude the singularities, one has: 



z,/l L J a —ioo v J 

xexp[~K(vcr {u)2 -c v a^)}, (Cll) 

where v = u + ao ■ Square completion and a change of 
variables w = 2v gives: 



Z {n \K) 



i/2 



K 



2ao+ioo 



dw 



2ao — ioo 



xexp[^nK(w + l/w)]w n l 2 

= ^K(27r/K) n / 2 I n/2 _ 1 (nK), (C12) 

which is an integral representation of the modified Bessel 
function of the first kind. Thus, the spin integrations 
can be done successively, the full partition function being 
given by: 

Z$>(K) = [(nK/2) 1 - n / 2 T(n/2)I n/2 _ 1 (nK)] N -\ 

where the T function enters through the normalization 
factor Affl and the last integration, which corresponds 
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to Z( n '(0). Considering the two-point function, an ad- 
ditional factor <Tj • <tj, i < j, is inserted in the inte- 
grand of Eq. (CI). Again starting the integration with 
the last spin cr^, the first N — j integrations are unal- 
tered. The integration over <jj gives addit ional factors 
of c„/2v from the Gaussian integration Eq. (C12), where 



now c u = Oj-ii so that one is left with 



K 



and the form of the integrand for the next integrations 
is unchanged. The integration over <Ti adds a factor of 



n since c v above becomes af ^ and J2 U U T > ' a 
followed by another i — 1 integrations of the partition- 
function type. With u = u(nK) = I n / 2 (nK) and v = 
v(nK) = I n /2~i{nK) one arrives at: 



» Jy) 



n. 



v N-l+i-j u j-i 



— (OTi ■ <Ta) = 



„N-1 



(u/v) 



3-i 



(C15) 



four-point function by analogy: 

\{{cn ■ * 3 ){cr k ■ *{)) = v N - l u l - k v k ^u j - i v i - 1 v 1 - N 
= (u/u) (l_fc)+(j_i) , (C16) 

where i < j < k < I is understood. For the special case 
of energy-energy correlations one has: 



1 



(Oi ■ tTi+l){tTj ■ O-j+i)) = (u/v) 2 



(C17) 



which does not depend on the distance \j — i\. Hence 
the connected energy-energy correlation function van- 
ishes exactly even for finite n in one dimension. The 
n — > oo limit of this expression is given by: 



1 4K 2 
— ((a-i ■ a i+ i)((Tj ■ (Tj+i)) = 



[1 + y/1 + (2K) 2 ] 2 ' 



(C18) 



From this it is straightforward to derive the form of the 



* Electronic address: weigel@itp.uni- leipzig.de 
t Electronic address: jankc Sitp.iini-lcipzig.de 

1 J. Zinn- Justin, Quantum Field Theory and Critical Phe- 
nomena (Clarendon Press, Oxford, 1996). 

2 L. P. Kadanoff, in Phase Transitions and Critical Phenom- 
ena, edited by C. Domb and M. S. Green (Academic Press, 
New York, 1976), Vol. 5A. 

3 M. E. Fisher and M. N. Barber, Phys. Rev. Lett. 28, 1516 
(1972). 

4 M. E. Fisher, Rev. Mod. Phys. 46, 597 (1974). 

5 M. E. Barber, in Phase Transitions and Critical Phenom- 
ena, edited by C. Domb and J. L. Lebowitz (Academic 
Press, New York, 1983), Vol. 8, p. 146. 

6 K. Binder, in Computational Methods in Field Theory, 
edited by H. Gausterer and C. B. Lang (Springer, Berlin, 
1992), p. 59. 

7 B. Diinweg, in Monte Carlo and Molecular Dynamics of 
Condensed Matter Systems, edited by K. Binder and C. 
Ciccotti (Italian Physical Society, Bologna, 1995). 

s A. A. Belavin, A. M. Polyakov, and A. B. Zamolodchikov, 

Nucl. Phys. B 241, 333 (1984). 
9 J. L. Cardy, Nucl. Phys. B 270, 186 (1986). 
10 J. L. Cardy, Nucl. Phys. B 275, 200 (1986). 

M. Henkel, Conformal Invariance and Critical Phenomena 
(Springer, Berlin/Heidelberg/New York, 1999). 
" J. L. Cardy, J. Phys. A 17, L385 (1984). 

R. B. Griffiths, Phys. Rev. Lett. 26, 1479 (1970). 
V. Privman and M. E. Fisher, Phys. Rev. B 30, 322 (1984). 
C. Domb, Adv. in Phys. 9, 1 (1960). 
M. P. Nightingale, Physica A 83, 561 (1976). 
B. Derrida and L. De Seze, J. Physique 43, 475 (1982). 
P. Nightingale and H. W. J. Blote, J. Phys. A 16, L657 
(1983). 

M. Weigel, Diploma thesis, Johannes-Gutenberg- 
Universitat Mainz, 1998, unpublished. 
J. L. Cardy, J. Phys. A 17, L961 (1984). 



21 J. L. Cardy, in Phase Transitions and Critical Phenomena, 
edited by C. Domb and J. L. Lebowitz (Academic Press, 
London, 1987), Vol. 11, p. 55. 

22 J. M. Luck, J. Phys. A 15, L169 (1982). 

23 K. A. Penson and M. Kolb, Phys. Rev. B 29, 2854 (1984). 

24 J. L. Cardy, J. Phys. A 18, L757 (1985). 

25 M. Weigel and W. Janke, Europhys. Lett., in print. 

26 M. Henkel, J. Phys. A 19, L247 (1986). 

27 M. Henkel, J. Phys. A 20, L769 (1987). 

28 M. Henkel, J. Phys. A 20, 3969 (1987). 

29 R. A. Weston, Phys. Lett. B 248, 340 (1990). 

30 M. Weigel and W. Janke, Phys. Rev. Lett. 82, 2318 (1999). 

31 R. Guida and J. Zinn- Justin, J. Phys. A 31, 8103 (1998). 

32 P. Butera and M. Comi, Phys. Rev. B 56, 8212 (1997). 

33 M. Campostrini, A. Pelissetto, P. Rossi, and E. Vicari, 
Phys. Rev. E 60, 3526 (1999). 

34 A. M. Ferrenberg and D. P. Landau, Phys. Rev. B 44, 5081 
(1991). 

35 H. W. J. Blote, E. Luijten, and J. R. Herringa, J. Phys. A 
28, 6289 (1995). 

36 H. W. J. Blote, L. N. Shchur, and A. L. Talapov, Int. J. 
Mod. Phys. C 10, 1137 (1999). 

37 M. Hasenbusch, K. Pinn, and S. Vinti, Phys. Rev. B 59, 
11471 (1999). 

38 U. Wolff, Phys. Rev. Lett. 62, 361 (1989). 

39 U. Wolff, Nucl. Phys. B 334, 581 (1990). 

40 W. Janke and K. Nather, Phys. Rev. B 48, 7419 (1993). 

41 H. Flyvbjerg and H. G. Petersen, J. Chem. Phys. 91, 461 
(1989). 

42 B. Efron, The Jackknife, the Bootstrap and Other Resam- 
pling Plans (Society for Industrial and Applied Mathemat- 
ics [SIAM], Philadelphia, 1982). 

43 B. A. Berg, Comput. Phys. Commun. 69, 7 (1992). 

44 A. L. Talapov and H. W. J. Blote, J. Phys. A 29, 5727 
(1996). 

45 A. M. Ferrenberg and R. H. Swendsen, Phys. Rev. Lett. 



18 



61, 2635 (1988). 

46 A. M. Ferrenberg and R. H. Swendsen, Phys. Rev. Lett. 
63, 1658(E) (1989). 

47 M. Campostrini, A. Pelissetto, P. Rossi, and E. Vicari, 
Phys. Rev. B 61, 5905 (2000). 

48 P. Butera, M. Comi, and A. J. Guttmann, Phys. Rev. B 
48, 13987 (1993). 

49 P. Butera and M. Comi, Phys. Rev. B 52, 6185 (1995). 

50 W. Janke, Phys. Lett. A 148, 306 (1990). 

51 H. G. Ballesteros, L. A. Fernandez, V. Martin-Mayor, and 
A. Munoz Sudupe, Phys. Lett. B 387, 125 (1996). 

52 M. Hasenbusch and T. Torok, J. Phys. A 32, 6361 (1999). 

53 M. Hasenbusch and S. Meyer, Phys. Lett. B 241, 238 
(1990). 

54 A. P. Gottlob and M. Hasenbusch, Physica A 201, 593 
(1993). 

55 J. Adler, C. Holm, and W. Janke, Physica A 201, 581 
(1993). 

56 P. Peczak, A. M. Ferrenberg, and D. P. Landau, Phys. Rev. 
B 43, 6087 (1991). 

57 K. Chen, A. M. Ferrenberg, and D. P. Landau, Phys. Rev. 
B 48, 3249 (1993). 

58 C. Holm and W. Janke, Phys. Rev. B 48, 936 (1993). 

59 C. Holm and W. Janke, Phys. Lett. A 173, 8 (1993). 

60 P. Reinicke, J. Phys. A 20, 4501 and 5325 (1987). 

61 W. Janke, in Computational Physics, edited by K. H. Hoff- 



mann and M. Schreiber (Springer, Berlin, 1996), p. 10. 

62 S. A. Antonenko and A. I. Sokolov, Phys. Rev. E 51, 1894 
(1995). 

63 H. Kleinert, Phys. Rev. D 60, 085001 (1999). 

64 M. A. Yurishchev, Phys. Rev. B 50, 13533 (1994). 

65 M. A. Yurishchev, Phys. Rev. E 55, 3915 (1997). 

66 H. E. Stanley, Phys. Rev. 176, 718 (1968). 

67 M. Henkel, J. Phys. A 21, L227 (1988). 

68 M. Henkel and R. A. Weston, J. Phys. A 25, L207 (1992). 

69 S. Allen and R. K. Pathria, J. Phys. A 26, 5173 (1993). 

70 E. Brezin, J. Physique 43, 15 (1982). 

71 T. H. Berlin and M. Kac, Phys. Rev. 86, 821 (1952). 

72 E. Helfand, Phys. Rev. 183, 562 (1969). 

73 M. Kac and C. J. Thompson, Physica Norvegica 5, 163 
(1971). 

74 H. W. Lewis and G. H. Wannier, Phys. Rev. 88, 682 (1952). 

75 C. C. Yan and G. H. Wannier, J. Math. Phys. 6, 1833 
(1965). 

76 H. E. Stanley, Phys. Rev. 179, 570 (1969). 

77 C. J. Thompson, in Phase Transitions and Critical Phe- 
nomena, edited by C. Domb and M. S. Green (Academic 
Press, New York, 1972), Vol. 2, p. 177. 

78 Considergig a closed chain is technically much more intri- 
cate, cp.Lj, but, of course, gives the same results in the 
thermodynamic limit N —* oo. 



